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ABSTRACT 

We numerically studied the explosion of a supernova caused by supersonic jets 
present in its center. The jets are assumed to be generated by a magneto-rotational 
mechanism when a stellar core collapses into a neutron star. We simulated the process 
of the jet propagation through the star, jet breakthrough, and the ejection of the 
supernova envelope by the lateral shocks generated during jet propagation. The 
end result of the interaction is a highly nonspherical supernova explosion with two 
high-velocity jets of material moving in polar directions, and a slower moving, oblate, 
highly distorted ejecta containing most of the supernova material. 

The jet-induced explosion is entirely due to the action of the jets on the surrounding 
star and does not depend on neutrino transport or re-acceleration of a stalled shock. 
The jet mechanism can explain the observed high polarization of Type Ib,c and 
Type II supernovae, pulsar kicks, very high velocity material observed in supernova 
remnants, indications that radioactive material was carried to the hydrogen-rich layers 
in SN1987A, and some others observations that are very difficult or impossible to 
explain by the neutrino energy deposition mechanism. The breakout of the jet from 
a compact, hydrogen-deficient core may account for the gamma-ray bursts and radio 
outburst associated with SN 1998bw/GRB980425. 



Subject headings: supernovae: general; individual (SN 1998bw) — gamma rays: bursts 
— pulsars: general — ISM: jets and outflows 



1. Introduction 

Recent observations of core collapse supernovae provide increasing evidence that the core 
collapse process is intrinsically asymmetric: 

1) The spectra of these supernovae are significantly polarized indicating asymmetric envelopes 
(Mendez et al. 1988; Hoflich 1991, 1995; Jeffrey 1991; Trammel et al. 1993, Tran et al. 1997). 
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The degree of polarization tend to vary inversely with the mass of the hydrogen envelope, being 
maximum for Type Ib/c events with no hydrogen (Wang et al. 1996; Wang, Wheeler &; Hoflich 
1999; Wheeler, Hoflich & Wang 1999). 

2) After the explosion, neutron stars are observed with high velocities, up to 1000 
km s-i (Strom et al. 1995). 

3) Observations of SN 1987A showed that radioactive material was brought to the hydrogen 
rich layers of the ejecta very quickly during the explosion (Lucy 1988; Sunyaev et al. 1987, Tueller 
et al. 1991). 

4) The remnant of the Cas A supernova shows rapidly moving oxygen-rich matter outside 
the nominal boundary of the remnant (Fesen & Gunderson, 1996) and evidence for two oppositely 
directed jets of high- velocity material (Fesen 1999; Reed, Hester, & Winkler 1999). 

5) High velocity "bullets" of matter have been observed in the Vela supernova remnant 
(Taylor et al. 1993) 

Understanding the mechanism of producing supernovae explosions by core collapse is a 
physics problem that has challenged researchers for decades (Hoyle & Fowler 1960; Colgate & 
White 1966). The current most sophisticated calculations based on the neutrino energy deposition 
mechanism are multidimensional and involve the convection of the newly formed neutron star. 
These, however, but have failed to produce robust explosions (Herant et al. 1994; Burrows, 
Hayes & Fryxell 1995; Janka &: Miiller 1996; Mezzacappa et al 1997; Lichtcnstadt, Khokhlov & 
Wheeler 1999). Even when successful, these models do not explain why SN 1998bw produced the 
strongest radio source ever associated with a supernova, probably requiring a relativistic blast 
wave (Kulkarni et al. 1998), or account for a probable link between SN 1998bw and the 7-ray 
burst GRB980425 observed in the same general location in the same general time frame (Galama 
et al. 1998). 

The discovery of pulsars led to early considerations of the role of rotating magnetized 
neutron stars in the explosion mechanism (LeBlanc and Wilson 1970; Ostriker k, Gunn 1971; 
Bisnovatyi-Kogan 1971). LeBlanc and Wilson studied the magneto-rotational core collapse of a 
7Mq star. They numerically solved the two-dimensional MHD equations coupled to the equation 
for neutrino transport. Their simulations showed the formation of two oppositely directed, 
high-density, supersonic jets of material emanating from the collapsed core. They estimated 
that at the surface located ~ 4 x 10^ cm from the center, the jet carried away ~ 10^^ g with 
~ 1 — 2 X 10^^ ergs in ~ 1 s. The magnetic field generated in this calculation was ~ 10-*^^ Gauss. 
Evidence now exists for strongly magnetized neutron stars, "magnetars" (Duncan &; Thompson 
1992; Kouveliotou et al. 1998). 

The LeBlanc- Wilson mechanism is extremely asymmetric and contains jets. Their calculations 
only followed the jet to a distance of ~ 10^ cm, whereas a stellar core has a radius of 10^*^ cm 
or more. The issues that arise are: how can this asymmetry propagate to much larger distances 
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inside the star? Can these jets induce asymmetry at distances comparable to the stellar radius, or 
even push through the entire star and exit? 

In this paper, we model the explosion of a core collapse supernova assuming that the 
LeBlanc- Wilson mechanism has operated in the center. Wc take a ISMq main-sequence star 
evolved to the point of the explosion (Stranicro, Chieffi &: Limongi 1999) and assume that the star 
has lost all of its hydrogen envelope before the explosion. The resulting 4.IM0 model of a helium 
star corresponds to the explosion of a Type lb or Ic supernova. The simulations show that the 
jets cause a very asymmetric explosion of the star. Most of the observations of asymmetries listed 
above can be explained by this process. 

2. Numerical Simulations 

Figure 1 presents a schematic of the setup of the computation. The computational domain 
is a cube of size L = 1.5 x 10^^ cm with a spherical helium star of radius i?star = 1-88 x lO^'^ cm 
and mass Mgtar — 4.1Mq placed in the center. The distribution of physical parameters inside 
the star is shown in Figure 2. The innermost part with mass Mcore — I-OMq and radius 
-Rcore = 3.82 X 10^ cm, consisting of Fe and Si, is assumed to have collapsed on a timescale much 
faster than the outer, lower-density material. It is removed and replaced by a point gravitational 
source with mass Mcore representing the newly formed neutron star. The remaining mass, from 
2± 1.6 to ~ 4.IM0, consists of an 0-Ne-Mg inner layer surrounded by the C-0 and He-envelopes. 
This structure is mapped onto the computational domain from i?core to .Rstar- 

At -Rcore f'lid the outer boundary of the computational domain, we impose an outflow 
boundary condition assuming zero pressure, velocity, and density gradients. At two polar 
locations where the jets are initiated at Rcore, we impose an inflow with velocity vj, density pj 
and pressure Pj. The jet parameters are chosen to represent the results of LeBlanc & Wilson 
(1970). At -Rcore) the jet density and pressure are the same as those of the background material, 
pj = 6.5 X 10^ g cm~^ and Pj = 1.0 x 10^^ ergs cm~^, respectively. The radii of the cylindrical jets 
entering the computational domain are approximately rj = 1.2 x 10^ cm. 

For the first 0.5 s, the jet velocity at -Rcorc is kept constant at Vj = 3.22 x 10^ cm s"^. 
This results in a mass flux rate of ~ 9.5 x 10^^ g s~^ with an energy deposition rate 
dE/dt = 5 X 10^'^ crgs/s for each jet. After 0.5 s, the velocity of the jets at i?corc was gradually 
decreased to zero at approximately 2 s. The total energy deposited by the jets is £'j ~ 9 x lO^'^ ergs 
and the total mass ejected is Mj ~ 2 x 10^^ grams or ~ O.IM©. These parameters are consistent 
within, but somewhat less than, those of the LeBlanc- Wilson model. The amount of material 
ejected is less than that which falls through the inner boundary during the jet operation, 
~ 4 X 10^^ g. This amounts to an implicit assumption that ^ 1/2 of the matter accreted is 
channeled back out into the jets. More accurate jet parameters can only be determined by 
self-consistently modeling the formation of the jets in the vicinity of a neutron star. 
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The stellar material was described by the time-dependent, compressible, Euler equations 
for inviscid flow with an ideal gas equation of state P = E{'y — 1) with constant 7 = 5/3. 
The Euler equations were integrated using an explicit, second-order accurate, Godunov type, 
adaptive- mesh-refinement, massively parallel, Fully-Threaded Tree (FTT) program, ALLA 
(Khokhlov 1998, Khokhlov Sz Chtchelkanova 1999). Euler fluxes were evaluated by solving a 
Riemann problem at cell interfaces. FTT discretization of the computational domain allowed the 
mesh to be refined or coarsened at the level of individual cells. Physical scales involved in the 
simulation range from the size of the computational domain (1.5 x 10^^ cm) to the jet diameter 
(~ 10^ cm) and span at least three orders of magnitude. Wc used a cartesian, nonuniformly 
refined FTT mesh with fine cells Amin — 3.7 X lO'^ cm near .Rcore to resolve the jets, and with 
cell size increasing towards the outer boundary of the computational domain where the cell size 
was Amax = 2.3 X 10^ cm. This mesh was fixed from initial time to 6 s of physical time. After 
that, the inner parts were coarsened near the center by a factor of four, and the central hole 
was eliminated. At this moment, the jets have exited the star and the details of the flow near 
i?corc do not affect the essential features of the explosion. In this first, demonstration calculation, 
we did not use the time-adaptive mesh refinement capability of ALLA. It will be used to follow 
shocks and mixing processes with higher resolution in future simulations. We computed the entire 
configuration including both jets and assuming no symmetries. The total number of computational 
cells used in the simulation was ~ 2 x 10^, whereas a uniform resolution Amin would have required 
- 7 X 10^° cells. 



3. Results and Discussion 

Figure 3 shows the propagation of the jet inside the star. As the jets move outwards, they 
remain coUimated and do not develop much internal structure. A bow shock forms at the head of 
the jet and spreads in all directions, roughly cylindrically around each jet. The sound crossing time 
r(r) = r/as{r) is shown as a function of stellar radius r in Figure 2, where as{r) is the sound speed 
at a given radius for the initial stellar model. It might be expected that if energy were released 
at the center of a star on a timescale much shorter that r(r), the effect of energy deposition at r 
would resemble that of a strong point explosion. In particular, the jet characteristic time Tj- ~ 1 s 
is much shorter than the sound crossing time of the star, r(i?star) ~ 10^ s (Figure 2). Nonetheless, 
these jets stay coUimated enough to reach the surface as strong jets. 

It is known that supersonic jets stay coUimated for a long distance. For example, Norman ct 
al. (1983) simulated supersonic jets with densities pj both less than and greater than a uniform 
background, pi,. Jets with Pj/pb > 1 developed a bow shock and little internal structure. Our jets 
resemble those with pj/pb > 1- The stellar matter is shocked by the bow shock, and then flows 
out and acts as a high-pressure confining medium by forming a cocoon around the jet. The sound 
crossing time of the dense 0-Ne-Mg mantle, r(i? ~ 10^ cm) ~ 10 s, is only ten times longer than 
Tj, and the jets are capable of penetrating this dense inner part of the star in ~ 2 s. By the time 
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the jets penetrate into the less dense C-0 and He layers, the inflow of material into the jets has 
been turned off. By this time, however, the jets have become long bullets of high-density material 
moving through the background low-density material almost ballistically. The higher pressures 
in these jets cause them to spread laterally. This spreading is limited by a secondary shock that 
forms around each jet between the jet and the material already shocked by the bow shock. The 
radius of the jets, ^ 3 x 10^ cm as they emerge from the star, is larger than the initial radius, 
~ 10^ cm, but it is still significantly less than the radius of the star. 

After about 5.9 s, the bow shock reaches the edge of the star and breaks through. Figure 4 
shows the subsequent evolution of the star after the breakthrough. By c± 20 s, most of the material 
in the jets has left the star propagates into the interstellar medium ballistically. We estimate the 
total mass in these two jets as Mj » O.OSMq and the total kinetic energy Ej 2.5 x 10^'^ ergs. 
The average velocity of the jet is about 25,000 km s~^. 

The laterally expanding bow shocks generated by the jets (Figure 3) move towards the equator 
where they collide with each other. The collision of the shocks first produces a regular reflection 
that then becomes a Mach reflection. The Mach stem moves outwards along the equatorial plane. 
The result is that the material in the equatorial plane is compressed and accelerated more than 
material in other directions (excluding the jet material). At t ~ 29 s, the Mach stem reaches 
the outer edge of the star, and the star begins to settle into the free expansion regime. The 
computation was terminated at ~ 35 s, before free expansion was attained. The stellar ejecta at 
this time is highly asymmetric. The density contour of 50 g cm~^, which is the average density 
of the ejecta at this time, forms an oblate configuration with the equator-to-polar velocity ratio 
2± 2/1. Complex shock and rarefaction interactions inside the expanding envelope will continue 
to change the distribution of the parameters inside the ejecta. Nonetheless, we expect that the 
resulting configuration will resemble an oblate ellipsoid with a very high degree of asymmetry, axis 
ratios > 2. 



4. Conclusions 

We have numerically studied the explosion of a supernova caused by supersonic jets generated 
in the center of the supernova as a result of the core collapse into a neutron star. We simulated 
the process of the jet propagation through the star, jet breakthrough, and the ejection of the 
supernova envelope by the lateral shocks generated during jet propagation. The end result of the 
interaction is a highly nonspherical supernova explosion with two high- velocity jets of material 
moving in polar directions ahead of an oblate, highly distorted ejecta containing most of the 
supernova material. Below we argue that such a model explains many of the observations that are 
difficult or impossible to explain by the neutrino deposition explosion mechanisms. 

We have assumed that the jets were generated by a magneto-rotational mechanism during 
core collapse and neutron star formation (LeBlanc &: Wilson 1970). That collimated jets could 
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be a common phenomenon in core collapse supernovae and be associated with 7-ray bursts was 
raised recently by Wang Sz Wheeler (1998). A different mechanism of jet generation involving 
neutrino radiation during collapse of a very massive star into a black hole has been recently 
discussed by MacFayden & Woosley (1999) in the context of a "failed" supernova, also to explain 
7-ray bursts. Low density relativistic jets may also be produced by the intense radiation of the 
newly born pulsar, as discussed by Blackman &; Yi (1998) and Yi et al. (1999). We found in our 
preliminary simulations (not presented here) that lower density and higher velocity jets than the 
one considered in this paper may produce similar hydrodynamical effects. 

The asymmetric explosion generated in this calculation provides ejection velocities that are 
comparable to those observed in supernovae. For this particular calculation, an energy of 2.5 x 10^° 
ergs is invested in the jet and the star of ~ 2.5M0 is ejected with kinetic energy of 6.5 x 10^*^ ergs 
and average velocity 3,000 — 4,000 km s~^. Increasing the jet opening angle, jet duration, or jet 
velocity would result in a more powerful explosion. The density and velocity profiles of the main 
ejecta (excluding jets) are oblate with equator to polar ratios greater than 2/1. This structure will 
produce significant polarization, of order 1% or more as observed in bare-core supernovae (Hoflich, 
Wheeler & Wang 1999). 

The two polar jets move outward from the star with a speed ~ 25, 000 km s"^, much greater 
than the ejecta itself. They may be detected in supernova remnants and might account for the 
evidence of jets in in Cas A (Fesen & Gunderson 1996; Reed, Hester & Winkler 1999). 

The composition of the jets must reflect the composition of the innermost parts of the star, 
and should contain heavy and intermediate-mass elements. During the explosion, the jets would 
bring heavy and intermediate mass elements into the outer layers. This will influence the spectral 
and polarization properties of a supernova. Here we considered a bare helium core, but if the core 
were inside a hydrogen envelope, the explosion would remain very inhomogeneous. Radioactive 
elements could potentially be carried into the hydrogen envelope. This could explain the early 
appearance of X-rays, as in SN1987A. It is plausible that a sufficiently powerful jet could even 
penetrate a hydrogen envelope. 

We assumed that the jets are identical which is not the general case. Any momentum 
imbalance might impart a kick to the neutron star. From momentum conservation, we estimate 
the required difference between the inflow velocities of the jets, Avj, be of the order of 



where Vns is the kick velocity, and we have taken the neutron mass Mjv5 = l.SM©, and the jet 
mass Mj = 10^^ g. Although the required jet asymmetry, = 1, to produce a 1,000 km/s kick 
may seem extreme, the parameters of jets selected for this calculations are mild. If the duration 
of the jets is increased by a factor of two, asymmetry of only 0.5 would be required. 




l,000km/s 




When the jets break through the stellar photosphere, a small amount of mass will be 
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accelerated through the density gradient to very high velocities. Our resolution was not enough 

and the code does not have a relativistic Riemann solver incorporated to make quantitative 
predictions; however, a small fraction of the material at the stellar surface was observed to move 
with a velocity of up to ~ 90,000 km s~^. This may, in principle, lead to the 7-ray burst and the 
radio outburst similar to those associated with SN 1998bw/GRB980425. 

The jet-induced explosion of a supernova computed in this paper is entirely due to the 
action of the jet on the surrounding star. The mechanism that determines the energy of such an 
explosion must be related to the shut-off of the accretion onto the neutron star by the lateral 
shocks that accelerate the material outwards. The explosion thus does not depend on neutrino 
transport or re-acceleration of the stalled shock. This work opens many issues that require further 
investigation. A study must be made of different input parameters, including properties of the jets 
and of the initial star, and the jet engine mechanisms must be studied as well. These studies are 
currently underway. 
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are grateful to Rob Duncan and Insu Yi for helpful discussions, and to Almadena Chtchelkanova 
for the development of massively parallel software used in the simulations. This research was 
supported in part by NSF Grant 95-28110, NASA Grant NAG 5-2888, NASA Grant LSTA-98-022 
and a grant from the Texas Advanced Research Program. The Laboratory for Computational 
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Program for support. 
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Figure Captions 

Fig. 1. — Schematic of the simulation. 

Fig. 2. — Initial Conditions. The distribution of physical parameters inside the innermost 5Mq of 
the 15M0 stellar model of Straniero et al. (1999). The Fe-Si inner part is assumed to collapse into 
a neutron star. The 0-Ne-Mg, C-0, and He layers are mapped onto the computational domain 
(see Section 2). 

Fig. 3. — Jet propagation inside the star. The frames show the density in the x-z plane passing 
through the center of computational domain. Time since the beginning of the simulation is given 
in the upper left corner of each frame. The sizes of two upper frames are Ax = 7.2 x 10^ cm and 
Az = 9.0 X 10^ cm. The sizes of the lower frames Ax = 3.6 x 10^° cm and Az = 4.5 x lO^'^ cm. 

Fig. 4. — Jet evolution after breakout. The frames show the density in the x-z plane passing through 
the center of computational domain. Time since the beginning of the simulation is given in the upper 
left corner of each frame. The sizes of frames are Ax = 6.1 x 10^° cm and Az = 1.125 x 10^^ cm. 
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Figure 2 
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